基于变指数分数阶全变差和整数阶全变差的图像恢复算法

结合变指数全变差(totalvariation, TV)和整数阶TV,提出一种变分图像恢复算法。该变分问题的能量泛函主要分为三个部分:变指数p(x)的分数阶TV正则化项、整数阶TV正则化项和数据保真项。该模型中的指数p(x)是与图像的梯度信息有关的函数。在理论上,由于分数阶导数和整数阶导数的结合,使得所提方法不仅能有效地去除图像噪音,保护图像的边界高频信息,还能更好地保留图像的纹理细节等中低频信息,同时可以极大地消除图像处理中产生的阶梯效应和散斑效应。在模型的求解上,利用变分法可以简单地将极小化泛函的优化问题转化为梯度下降流方程。最后,通过模拟数据和真实数据对本文所提方法进行了验证。试验结果表明,该方法可以去除噪声的同时,有效保持边界和纹理细节,并且对噪声是鲁棒的,具有一定的实际应用价值。     

基于GPU的多维树搜索算法在管道流体仿真中的应用

为提高邻近粒子搜索速度,使水力压裂过程中注入管道内的水流现象更好地满足实时性要求,在图形处理器(GPU)上基于光滑粒子流体动力学(SPH)算法,将结合GPU的并行操作能力和多维树搜索算法优点的多维树搜索算法应用在水力压裂管道流体模拟中。首先,给出了流体粒子的受力分析及更新颗粒位置的方法;然后,给出了该算法在管道流体模拟中的具体实现过程。最后,通过实验近似模拟水力压裂过程中管道与水流交互过程,并与其他算法进行对比分析。结果表明将本文算法应用于大规模管道流体模拟中,执行速度快了几倍,实现了水流模拟在管道中的实时性。     

箭型排布的矩形微槽平板上液膜流动特性

光滑平板降膜受表面张力和接触角的影响易收缩成溪流,导致传热表面出现干斑,为解决这一问题,提出箭型排布的矩形微槽平板。通过可靠的computational fluid dynamics（CFD）计算模拟两相流理论,建立三维非稳态平板降膜数学模型,研究了箭型排布的矩形微槽平板上的液膜流动特性,并探究了微槽宽度、深度及箭型夹角对液膜在平板上铺展效果的影响。结果表明：箭型排布的矩形微槽可有效增大液膜在平板横向的铺展面积,使液膜润湿面积增大,减少平板表面干斑;在120°箭型结构下,矩形微槽最优参数为宽0.5 mm,深0.3 mm,此时可将比湿面积由光板表面的62%~89%提高到84%~94%;低雷诺数时,120°箭型结构对液膜横向铺展引导效果显著,雷诺数增大时,90°箭型结构引导效果更好。     

列车绕流的瞬态与稳态数值模拟对比

采用大涡模拟(Large Eddy Simulation,LES)和基于雷诺平均(Reynolds Average Navier-Stocks,RANS)的SST(Shear Stress Transport)k-ω湍流模型,分别对高速列车单车明线运行进行瞬态和稳态的仿真计算,通过与实车测试数据比较对数值模拟进行了验证.对比分析LES和RANS的计算结果发现:对于车头表面测压点,LES和RANS都能给出高精度的计算结果,且LES的瞬态计算结果表明,表面压力最大值在一个很宽的范围内波动;对于列车绕流结构,LES较RANS表现出更强的小尺度涡的捕捉能力,尤其表现在复杂的尾流区;通过气动力系数的傅里叶变换分析了波动的频域特性.LES在较复杂列车模型外流场模拟中的高计算精度,及其广泛的结果信息可以为列车的系统耦合设计提供可靠的数据参考.     

基于灵敏度阻抗矩阵修正法的分层前推回代潮流算法

分布式电源(DG)接入配电网对潮流产生重要影响.本文分析了潮流计算中各种分布式电源模型及处理方法,引入灵敏度阻抗矩阵修正法更新PV节点的注入无功功率,结合辐射型配电网的特点,提出一种基于灵敏度阻抗矩阵修正法的分层前推回代潮流算法.该算法解决了前推回代潮流算法处理PV节点失效的问题,同时适用于含各类分布式电源的潮流计算.最后对含各种类型分布式电源的IEEE 33节点配电网进行潮流计算仿真,仿真结果验证了提出算法的有效性和快速性,并通过不同算例验证了算法的稳定性.     

铀矿山下垫面粗糙度对大气风场影响的数值模拟研究

以某实际铀矿山为对象,运用CFD方法,基于N-S方程建立三维大气风场的物理模型,研究三种下垫面粗糙度下（z=0.1 m,z=0.5 m,z=1.0 m）的大气风场结构.结果显示,随着粗糙度增大,近地面风场速度变化达到19%;随着高度上升,大约在200 m后,粗糙度影响减弱,风场等值线趋于平稳;风场流速越大,近壁面处粗糙度对风场的阻碍作用越明显.在粗糙度较大（z≥0.5 m）情况下,可适当通过提高污染物排放高度来加强大气的输运能力,以缓解局部地区的污染。     

基于Logistic回归模型的IT企业财务危机预警的实证研究

以我国信息技术上市的ST公司和非ST公司作为研究样本,运用因子分析以及Logistic回归分析的方法,综合传统的财务指标和现金流信息指标构建适合IT企业的财务危机预测模型.并对模型的预测效果进行回判,验证了该模型对IT行业具有很好的预警效果.     

两种结构型式控流筛管内流动特性的数值研究

控流筛管可通过控制油井部分井段产出液流量实现水平段同时出水。以计算流体动力学为基础,分别对导流式控件和迷宫式控件流道内流体的流动特征进行数值模拟,并对两种控件的适用范围进行对比分析。数值计算结果表明,两种控件由流道入口到出口的流道压力均呈线性降低;导流式控件流道中的流动近似呈平推流,涡流和湍动不显著;低黏度条件下,迷宫式控件流道内存在明显的一次、二次涡流,过流阻力显著增加,控流效果显著。过流流体黏度为25~1000 mPa·s时,导流式控件的性能较优;其他黏度范围内,迷宫式控件的性能较优。     

基于模拟有限差分法的离散裂缝模型两相流动模拟

模拟有限差分作为一种新型数值计算方法,因其良好的局部守恒性和对复杂网格系统的适用性,在计算流体力学和油藏数值模拟中得到应用。将模拟有限差分方法进一步应用于离散裂缝流动模拟中,对模拟有限差分法的基本原理进行详细阐述,建立相应的离散裂缝数值计算格式,并采用IMPES方法对其两相流问题进行求解,并与实验结果对比。对比结果验证了新方法和程序的正确性,通过复杂离散裂缝算例进一步验证了方法的正确性和程序的鲁棒性。     

局部流水FFT处理器设计

讨论局部流水FFT处理器中的两个主要模块:蝶形运算流水线和地址产生器的设计.基于对基2蝶形单元的"深"反馈,提出一种称之为R2SD2 F(radix-2single"deep"delay feedback,基2单路深度延时反馈)的流水线结构.该流水线中的蝶形处理单元仅由两个复数加法器组成,可以工作在基4/基2/直通三种模式下,因此由两个如此蝶形处理单元组成的R2SD2F流水线可以在一次循环中选择完成基16/基8/基4/基2运算.在完成长为N(假定N为4的整数次幂)点的DFT运算时,该流水线所需的主要硬件有log4N-1个复数乘法器和2log4N个复数加法器.作为一个整体,给出局部流水FFT处理器中的地址产生方法和旋转因子存取结构.